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The magnetization dynamics of ferromagnets are often formulated in terms of the Landau-Lifshitz- 
Gilbert (LLG) equation. The reactive part of this equation describes the response of the magnetiza- 
tion in terms of effective fields, whereas the dissipative part is parameterized by the Gilbert damping 
tensor. We formulate a scattering theory for the magnetization dynamics and map this description 
on the linearized LLG equation by attaching electric contacts to the ferromagnet. The reactive part 
can then be expressed in terms of the static scattering matrix. The dissipative contribution to the 
low-frequency magnetization dynamics can be described as an adiabatic energy pumping process 
to the electronic subsystem by the time-dependent magnetization. The Gilbert damping tensor 
depends on the time derivative of the scattering matrix as a function of the magnetization direction. 
By the fluctuation-dissipation theorem, the fluctuations of the effective fields can also be formulated 
in terms of the quasistatic scattering matrix. The theory is formulated for general magnetization 
textures and worked out for monodomain precessions and domain wall motions. We prove that the 
Gilbert damping from scattering theory is identical to the result obtained by the Kubo formalism. 

PACS numbers; 75.40.Gb,76.60.Es,72.25.Mk 



o 



> 

O 



X 



I. INTRODUCTION 

Ferromagnets develop a spontaneous magnetization 
below the Curie temperature. The long-vifavelength mod- 
ulations of the magnetization direction consist of spin 
waves, the low-lying elementary excitations (Goldstonc 
modes) of the ordered state. When the thermal energy is 
much smaller than the microscopic exchange energy, the 
magnetization dynamics can be phenomenologically ex- 
pressed in a generalized Landau-Lifshitz-Gilbert (LLG) 
form: 

m(r, t) = -7m(r, t) x [Heff (r, t) + h(r, t)] + 

m{r,t)x y"dr' [a [m] (r,r')m(r',i)] , (1) 

where the magnetization texture is described by m(r, t), 
the unit vector along the magnetization direction at po- 
sition r and time t, m{r,t) ~ dm.{r, t)/dt, 7 = g^s/fi is 
the gyromagnetic ratio in terms of the (7-factor (sa 2 for 
free electrons) and the Bohr magneton ^b- The Gilbert 
damping a is a nonlocal symmetric 3x3 tensor that is 
a functional of m. The Gilbert damping tensor is com- 
monly approximated to be diagonal and isotropic (i), lo- 
cal (1), and independent of the magnetization m, with 
diagonal elements 



aii(r, r') = aS{r — r'). 



(2) 



The linearized version of the LLG equation for small- 
amplitude excitations has been derived microscopicallyii 



It has been used very successfully to describe the mea- 
sured response of ferromagnetic bulk materials and thin 
films in terms of a small number of adjustable, material- 
specific parameters. The experiment of choice is fer- 
romagnetic resonance (FMR), which probes the small- 
amplitude coherent precession of the magnet 1^ The 
Gilbert damping model in the local and time-independent 
approximation has important ramifications, such as a lin- 
ear dependence of the FMR line width on resonance fre- 
quency, that have been frequently found to be correct. 
The damping constant is technologically important since 
it governs the switching rate of ferromagnets driven by 
external magnetic fields or electric currentsj^ In spatially 
dependent magnetization textures, the nonlocal charac- 
ter of the damping can be significant as well4r— Moti- 
vated by the belief that the Gilbert damping constant is 
an important material property, we set out here to under- 
stand its physical origins from first principles. We focus 
on the well studied and technologically important itiner- 
ant ferromagnets, although the formalism can be used in 
principle for any magnetic system. 

The reactive dynamics within the LLG Eq. ([TJ is de- 
scribed by the thermodynamic potential f2[M] as a func- 
tional of the magnetization. The efi'ective magnetic field 
Hoff[M](r) = — (517/(5M(r) is the functional derivative 
with respect to the local magnetization M(r) = Msm(r), 
including the external magnetic field Hoxt, the magnetic 
dipolar field Hd, the texture-dependent exchange energy, 
and crystal field anisotropics. Mg is the saturation mag- 
netization density. Thermal fluctuations can be included 
by a stochastic magnetic fleld h(r, t) with zero time av- 
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FIG. 1: Schematic picture of a ferromagnet (F) in contact 
with a thermal bath (reservoirs) via metaUic normal metal 
leads (N). 

erage, (h) = 0, and white-noise correlation>i 

{h,ir,t)h,ir',t')) = ^a,, [m] (r,r')<5(t - t'), (3) 

where Ms is the magnetization, i and j are the Cartesian 
indices, and T is the temperature. This relation is a con- 
sequence of the fluctuation-dissipation theorem (FDT) in 
the classical (Maxwell-Boltzmann) limit. 

The scattering {S-) matrix is defined in the space of 
the transport channels that connect a scattering region 
(the sample) to real or fictitious thermodynamic (left 
and right) reservoirs by electric contacts with leads that 
are modeled as ideal wave guides. Scattering matri- 
ces are known to describe transport properties, such as 
the giant magnetoresistance, spin pumping, and current- 
induced magnetization dynamics in layered normal-metal 
(N) I ferromagnet (F)<^"— When the ferromagnet is part 
of an open system as in Fig. [1] also f2 can be expressed 
in terms of the scattering matrix, which has been used 
to express the non-local exchange coupling between fer- 
romagnetic layers through conducting spacersJ^ We will 
show here that the scattering matrix description of the 
effective magnetic fields is valid even when the system is 
closed, provided the dominant contribution comes from 
the electronic band structure, scattering potential disor- 
der, and spin-orbit interaction. 

Scattering theory can also be used to compute the 
Gilbert damping tensor a for magnetization dynamicsJ^ 
The energy loss rate of the scattering region can be ex- 
pressed in terms of the time-dependent S'-matrix. To this 
end, the theory of adiabatic quantum pumping has to be 
generalized to describe dissipation in a metallic ferromag- 
net. The Gilbert damping tensor is found by evaluating 
the energy pumping out of the ferromagnet and relat- 
ing it to the energy loss that is dictated by the LLG 
equation. In this way, it is proven that the Gilbert phe- 
nomenology is valid beyond the linear response regime 
of small magnetization amplitudes. The key approxima- 
tion that is necessary to derive Eq. ([T]) including a is the 
(adiabatic) assumption that the ferromagnetic resonance 
frequency ojfmr that characterizes the magnetization dy- 
namics is small compared to internal energy scale set by 
the exchange splitting A and spin-flip relaxation rates 
Tg. The LLG phenomenology works well for ferromag- 
nets for which wfmr A/fi,, which is certainly the case 
for transition metal ferromagnets such as Fe and Co. 

Gilbert damping in transition-metal ferromagnets is 
generally believed to stem from the transfer of energy 
from the magnetic order parameter to the itinerant quasi- 



particle continuum. This requires either magnetic disor- 
der or spin-orbit interactions in combination with impu- 
rity/phonon scatteringi^ Since the heat capacitance of 
the ferromagnet is dominated by the lattice, the energy 
transferred to the quasiparticles will be dissipated to the 
lattice as heat. Here we focus on the limit in which clas- 
tic scattering dominates, such that the details of the heat 
transfer to the lattice docs not affect our results. Our ap- 
proach formally breaks down in sufficiently clean samples 
at high temperatures in which inelastic electron-phonon 
scattering dominates. Nevertheless, quantitative insight 
can be gained by our method even in that limit by mod- 
elling phonons by frozen deformations 

In the present formulation, the heat generated by the 
magnetization dynamics can escape only via the contacts 
to the electronic reservoirs. By computing this heat cur- 
rent through the contacts we access the total dissipa- 
tion rate. Part of the heat and spin current that es- 
capes the sample is due to spin pumping that causes 
energy and momentum loss even for otherwise dissipa- 
tion less magnetization dynamics. This process is now 
well understoodji^ For sufficiently large samples, the spin 
pumping contribution is overwhelmed by the dissipation 
in the bulk of the ferromagnet. Both contributions can 
be separated by studying the heat generation as a func- 
tion of the length of a wire. In principle, a voltage can be 
added to s tudy d issipation in the presence of electric cur- 
rents as in ll3lll4l . but we concentrate here on a common 
and constant chemical potential in both reservoirs. 

Although it is not a necessity, results can be simpli- 
fied by expanding the S'-matrix to lowest order in the 
amplitude of the magnetization dynamics. In this limit 
scattering theory and the Kubo linear response formal- 
ism for the dissipation can be directly compared. We 
will demonstrate explicitly that both approaches lead to 
identical results, which increases our confidence in our 
method. The coupling to the reservoirs of large samples 
is identified to play the same role as the infinitesimals in 
the Kubo approach that guarantee causality. 

Our formalism was introduced first in Ref. [ill lim- 
ited to the macrospin model and zero temperature. An 
extension to the friction associated with domain wall mo- 
tion was given in Ref. [l^ Here we show how to handle 
general magnetization textures and finite temperatures. 
Furthermore, we offer an alternative route to derive the 
Gilbert damping in terms of the scattering matrix from 
the thermal fluctuations of the effective field. We also 
explain in more detail the relation of the present theory 
to spin and charge pumping by magnetization textures. 

Our paper is organized in the following way. In Sec- 
tion HH wc introduce our microscopic model for the fer- 
romagnet. In Section IIIIl dissipation in the Landau- 
Lifshitz-Gilbert equation is exposed. The scattering the- 
ory of magnetization dynamics is developed in Sec. IIVI 
We discuss the Kubo formalism for the time-dependent 
magnetizations in Scc.|Vl before concluding our article in 
Sec. lVIl The Appendices provide technical derivations of 
spin, charge, and energy pumping in terms of the scat- 
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tcring matrix of the system. 

II. MODEL 

Our approach rests on density-functional theory 
(DFT), which is widely and successfully used to describe 
the electronic structure and magnetism in many fer- 
romagncts, including transition-metal ferromagnets and 
ferromagnetic semiconductors.— In the Kohn-Sham im- 
plementation of DFT, noninteracting hypothetical par- 
ticles experience an effective exchange-correlation poten- 
tial that leads to the same ground-state density as the in- 
teracting many-electron system.— A simple yet successful 
scheme is the local-density approximation to the effective 
potential. DFT theory can also handle time-dependent 
phenomena. We adopt here the adiabatic local-density 
approximation ( ALDA) , i. e. an exchange-correlation po- 
tential that is time-dependent, but local in time and 
spaceJ^ii^ As the name expresses, the ALDA is valid 
when the parametric time-dependence of the problem is 
adiabatic with respect to the electron time constants. 
Here we consider a magnetization direction that varies 
slowly in both space and time. The ALDA should be 
suited to treat magnetization dynamics, since the typical 
time scale (ipMR ~ 1/ (10 GHz) ^ 10~^°s) is long com- 
pared to the that associated with the Fermi and exchange 
energies, 1 — 10 eV leading to /i/A ~ 10^^'^s in transition 
metal ferromagnets. 

In the ALDA, the system is described by the time- 
dependent effective Schrodinger equation 

^ALDA*(r,<) =in|-*(r,<), (4) 
at 

where ^'(r,t) is the quasiparticle wave function at posi- 
tion r and time t. We consider a generic mean- field elec- 
tronic Hamiltonian that depends on the magnetization 
direction -ffALDA [m] and includes the periodic Hartree, 
exchange and correlation potentials and relativistic cor- 
rections such as the spin-orbit interaction. Impurity scat- 
tering including magnetic disorder is also represented by 
Halba- The magnetization m is allowed to vary in time 
and space. The total Hamiltonian depends additionally 
on the Zeeman energy of the magnetization in external 
Hcxt and dipolar magnetic fields: 

H = Halda [m] -Ah j drm- (He.t + H^) . (5) 

For this general Hamiltonian ([5]), our task is to de- 
duce an expression for the Gilbert damping tensor a. To 
this end, from the form of the Landau-Lifshitz-Gilbert 
equation ([3]) , it is clear that we should seek an expansion 



in terms of the slow variations of the magnetizations in 
time. Such an expansion is valid provided the adiabatic 
magnetization precession frequency is much less than the 
exchange splitting A or the spin-orbit energy which gov- 
erns spin relaxation of electrons. We discuss first dissi- 
pation in the LLG equation and subsequently compare 
it with the expressions from scattering theory of electron 
transport. This leads to a recipe to describe dissipation 
by first principles. Finally, we discuss the connection to 
the Kubo linear response formalism and prove that the 
two formulations are identical in linear response. 



III. DISSIPATION AND 
LANDAU-LIFSHITZ-GILBERT EQUATION 

The energy dissipation can be obtained from the solu- 
tion of the LLG Eq. ^ as 

E = -Af.y" dr[m(r,<)-Heff(r,t)] (6) 
= J '^^ J '^r'rii(r) • a [m] (r, r') • rh(r'). (7) 

The scattering theory of magnetization dissipation can be 
formulated for arbitrary spatiotemporal magnetization 
textures. Much insight can be gained for certain special 
cases. In small particles or high magnetic fields the col- 
lective magnetization motion is approximately constant 
in space and the "macrospin" model is valid in which 
all spatial dependences are disregarded. We will also 
consider special magnetization textures with a dynamics 
characterized by a number of dynamic (soft) collective 
coordinates ^a{t) counted by aiMi^ 

m(r,t) = m3t(r;Ua(i)}), (8) 

where nist is the profile at t — > —oo. This representation 
has proven to be very effective in handling magnetiza- 
tion dynamics of domain walls in ferromagnetic wires. 
The description is approximate, but (for few variables) 
it becomes exact in special limits, such as a transverse 
domain wall in wires below the Walker breakdown (see 
below); it becomes arbitrarily accurate by increasing the 
number of collective variables. The energy dissipation to 
lowest (quadratic) order in the rate of change of the 
collective coordinates is 

E=-J2^abiaL (9) 
ab 

The (symmetric) dissipation tensor Tab reads 
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f - ^ / dr 

1 



dr — „ , a [mj (r, r ) 



(10) 



r 



The equation of motion of the coUective coordinates un- 
der a force 
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(11) 



(12) 



introducing the antisymmetric and time-independent gy- 
rotropic tensor: 



■qab 
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dv mst(r) 



5mst(r) 5mst(r) 



96 



(13) 



We show below that 5^ and f can be expressed in terms 
of the scattering matrix. For our subsequent discussions 
it is necessary to include a fluctuating force f(t) (with 
{fit)) = 0), which has not been considered in Refs . [2011211 . 
From Eq. (jSj if follows the time correlation of f is white 
and obeys the fluctuation-dissipation theorem: 



{Ut)W))^2kBTfab5{t-t'). 



(14) 



In the following we illustrate the collective coordinate 
description of magnetization textures for the macrospin 
model and the Walker model for a transverse domain 
wall. The treatment is easily extended to other rigid 
textures such as magnetic vortices. 



A. Macrospin excitations 

When high magnetic flelds are applied or when the 
system dimensions are small the exchange stiffness dom- 
inates. In both limits the magnetization direction and 
its low energy excitations lie on the unit sphere and its 
magnetization dynamics is described by the polar angles 
e{t) and ip{t): 



m = (sin cos Lp, sin 9 sin cos ( 



(15) 



The diagonal components of the gyrotropic tensor vanish 
by (anti)symmetry fjee = 0, fj^^ = 0. Its off-diagonal 
components are 

= sm = -ilu,e. (16) 

7 

V is the particle volume and M^V the total magnetic 
moment. We now have two coupled equations of motion 



MsV 

7 

MsV ^ 
7 



'do 

an 



^Sine- — - {Tgg9 + Tg^^] =0, (17) 



sine* - — - {r^eO + T^ 



dip 



ipip'-P 



= 0. 



The thermodynamic potential D, determines the ballistic 
trajectories of the magnetization. The Gilbert damping 
tensor Tab will be computed below, but when isotropic 
and local, 



f = I5{r - r')Msa/j, 



(18) 



where 1 is a unit matrix in the Cartesian basis and a 
is the dimensionless Gilbert constant, Tgg ~ MgVa/j, 



r, 



0^p 







lp0 



and r,. 



'0MsVa/"f. 



B. Domain Wall Motion 

We focus on a one-dimensional model, in which the 
magnetization gradient, magnetic easy axis, and external 
magnetic field point along the wire (z) axis. The mag- 
netic energy of such a wire with transverse cross section 
S can be written as^ 



(19) 



n = AIsS / dz(j}{z) 



in terms of the one-dimensional energy density 



A 
2" 



dm 



dz 



Ham, + ^{l-mi) + ^mi, (20) 



where Ha is the applied field and A is the exchange stiff- 
ness. Here the easy-axis anisotropy is parametrized by 
an anisotropy constant Ki. In the case of a thin film 
wire, there is also a smaller anisotropy energy associated 
with the magnetization transverse to the wire governed 
by if 2- In a cylindrical wire from a material without 
crystal anisotropy (such as permalloy) K2 = 0. 

When the shape of such a domain wall is pre- 
served in the dynamics, three collective coordinates 
characterize the magnetization texture: the domain 
wall position ^i{t) = r^{t), the polar angle ^2(0 = 
(Pw{t), and the domain wall width Aw(i)- We con- 
sider a head-to-head transverse domain wall (a tail- 
to-tail wall can be treated analogously). m(z) = 
(sin 6'w cos ipvt , sin 6^ sin , cos 6'w) , where 



cos 0w = tanh ■ 



and 



esc 6'w = cosh ■ 



A. 



(21) 



(22) 



minimizes the energy ([^0]) under the constraint that the 
magnetization to the far left and right points towards the 
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domain wall. The off-diagonal elements are then fjri = 
= fjij. and fir^ = —2AIs/j = —fj^pj.- The energy ([^0]) 
reduces to 

n = M,S [A/X„ - 2Har + KiX^ + K2K cos^ <p^] . 

(23) 

Disregarding fluctuations, the equation of motion Eq. 
p2)) can be expanded as: 

2rw + a^^tp + a^r'f'w + a^xX^, = lK2^w sin 2tp^, 

(24) 

-2ip + arr-Tw + Oir^if + arxXv, ^ i'yHa, (25) 

A/A^ + axrr^ + ax^p + q^aaAw = + A'2 cos^ (^w, 

(26) 

where Ofah = ^Tab/MsS. 

When the Gilbert damping tensor is isotropic and local 
in the basis of the Cartesian coordinates, F = 15 {r — 
T')Msah 

2a tP'ol 
otrr = T— ; OL^p^p = 2a.X^; ax\ = (27) 
Aw dAw 

whereas all off-diagonal elements vanish. 

Most experiments are carried out on thin film ferro- 
magnetic wires for which K2 is finite. Dissipation is es- 
pecially simple below the Walker threshold, the regime 
in which the wall moves with a constant drift velocity, 
<^w = and^^ 



-2-^11 a I OLr 



(28) 



The Gilbert damping coefficient a^r can be obtained di- 
rectly from the scattering matrix by the parametric de- 
pendence of the scattering matrix on the center coordi- 
nate position When the Gilbert damping tensor is 
isotropic and local, we find — X^^Haj a. The domain 
wall width Aw = ^ Aj (K\ + Ki cos^ (/^w) and the out- 
of-plane angle (/3w = 5 arcsin27i/a/a^2- At the Walker- 
breakdown field {Ha)^^^ = aK2/ (27) the sliding domain 
wall becomes unstable. 

In a cylindrical wire without anisotropy, K2 — 0, is 
time-dependent and satisfies 



(2 -f a^r) . 



IV. SCATTERING THEORY OF MESOSCOPIC 
MAGNETIZATION DYNAMICS 

Scattering theory of transport phcnomena^i has 
proven its worth in the context of magnetoelectronics. 
It has been used advantageously to evaluate the non- 
local exchange interactions multilayers or spin valves,— 
the giant magnetoresistance?^ spin-transfer torque?^ and 
spin pumpingjiS We first review the scattering theory 
of equilibrium magnetic properties and anisotropy fields 
and then will turn to non-equilibrium transport. 



A. Conservative forces 

Considering only the electronic degrees of freedom in 
our model, the thermodynamic (grand) potential is de- 
fined as 



n = -fcRTluTre-^^*^"*-''^) 



(32) 



while /i is the chemical potential, and N is the number 
operator. The conservative force 



(33) 



can be computed for an open systems by defining a scat- 
tering region that is connected by ideal leads to reservoirs 
at common equilibrium. For a two-terminal device, the 
flow of charge, spin, and energy between the reservoirs 
can then be described in terms of the S'-matrix: 



r t' 
t r' 



(34) 



(29) 



where r is the matrix of probability amplitudes of states 
impinging from and reflected into the left reservoir, while 
t denotes the probability amplitudes of states incoming 
from the left and transmitted to the right. Similarly, 
r' and i' describes the probability amplitudes for states 
that originate from the right reservoir, r, r', i, and i' are 
matrices in the space spanned by eigenstates in the leads. 
We are interested in the free magnetic energy modulation 
by the magnetic configuration that allows evaluation of 
the forces Eq. (|33p . The free energy change reads 



while 



V "fv J 



ar 



(30) 



For isotropic and local Gilbert damping coefficients;^ 



1 + q;2 ' 



(31) 



In the next section, we formulate how the Gilbert scatter- 
ing tensor can be computed from time-dependent scat- 
tering theory. 



An = -keT / deAn{e) In I 



(35) 



where An(e)de is the change in the number of states at 
energy e and interval de, which can be expressed in terms 
of the scattering matriji^ 



An(e) 



1 d 



TrlnS'(e) 



27ri de 

Carrying out the derivative, we arrive at the force 



(36) 



(37) 
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where /(e) is the Fermi-Dirac distribution function with 
chemical potential ^. This estabhshed result will be re- 
produced and generalized to the description of dissipation 
and fluctuations below. 



B. Gilbert damping as energy pumping 



In the macrospin model the Gilbert damping tensor can 
then be expressed as 



de 



dl 
de 



Tr 



dSje) dS^e) 
drrii druj 



(41) 



Here we interpret Gilbert damping as an energy pump- 
ing process by equating the results for energy dissipa- 
tion from the microscopic adiabatic pumping formalism 
with the LLG phenomenology in terms of collective co- 
ordinates, Eq. The adiabatic energy loss rate of a 
scattering region in terms of scattering matrix at zero 
temperature has been derived in Refs. [26ll27l . In the ap- 
pendices, we generalize this result to finite temperatures: 



dt 



dt 



(38) 



Since we employ the adiabatic approximation, S'(e, t) is 
the energy-dependent scattering matrix for an instanta- 
neous ("frozen") scattering potential at time t. In a mag- 
netic system, the time dependence arises from its magne- 
tization dynamics, S{e,t) = S[m{t)]{e). In terms of the 
collective coordinates ^(<), S(e,t) = 5(e, {^(t)}) 



dS[m{t)] 
dt 



(39) 



where the approximate sign has been discussed in the 
previous section. We can now identify the dissipation 
tensor PH)) in terms of the scattering matrix 



Tab = 



47r 



dt - 



df 



Tr 



dSje) dS^e 



(40) 



where nii is a Cartesian component of the magnetization 
direction.. 



Gilbert damping and fluctuation-dissipation 
theorem 



At finite temperatures the forces acting on the mag- 
netization contain thermal fluctuations that are related 
to the Gilbert dissipation by the fluctuation-dissipation 
theorem, Eq. ([Ti]) . The dissipation tensor is therefore ac- 
cessible via the stochastic forces in thermal equilibrium. 

The time dependence of the force operators 



m 



(42) 



is caused by the thermal fluctuations of the magneti- 
zation. It is convenient to rearrange the Hamiltonian 
Halba into an unperturbed part that does not de- 
pend on the magnetization and a scattering potential 
^?ALDA(m) = ^^0 + V{m.). In the basis of scattering 
wave functions of the leads, the force operator reads 



I de I de'{ea\^\e'(i)di{e)dp{e')e^^^-^'^'/\ 



(43) 



r 



where djs annihilates an electron incident on the scatter- 
ing region, /3 labels the lead (left or right) and quantum 
numbers of the wave guide mode, and |e'/3) is an associ- 
ated scattering eigenstate at energy e'. We take again the 
left and right reservoirs to be in thermal equilibrium with 
the same chemical potentials, such that the expectation 
values 



{diie)d,{e'))^S^,S{e^e')fie). 



(44) 



The relation between the matrix element of the scattering 
potential and the S'-matrix 



S\e) 



dS{e) 



dV 

= -2■K^{ea\ — \eP) 

al3 '^^ 



(45) 



follows from the relation derived in Eq. ([BT]) below as 
well as unitarity of the 5- matrix, S = 1. Taking these 
relations into account, the expectation value oi'S is found 
to be Eq. ([37]) • We now consider the fluctuations in the 
force f{t) = ^{t) — {^{t)), which involves expectation 
values 

(4i (ei)a/Ji (e'i)aL (e2)a& (4)) 

- (ali(ei)a/3i(ei)) (oL^ (£2)0/32 (4)) 

= Sa^fjJ (ei - e'2) 5(3^a^5 (e'l - £3) /(ci) [1 - f{^2)] , 

(46) 

where we invoked Wick's theorem. Putting everything 
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together, we finaUy find 

{fa{t)fb{t')) ^2kBTS{t-t')Tab, 



(47) 



where Tab has been defined m Eq. (|40p . Comparing witli 
Eq. ((HI), we conclude that the dissipation tensor Tab 
governing the fluctuations is identical to the one obtained 
from the energy pumping, Eq. ()40|) . thereby confirming 
the fluctuation-dissipation theorem. 



consists of an equilibrium contribution and a term linear 
in the perturbed magnetization direction: 



datl) (t) = (daH)^^ +Y,j dt'Xabit - t')SUt') ■ 

(52) 

Here, we introduced the retarded susceptibility 



KUBO FORMULA 



Xabit - t') = --9{t - t') ^ [daH{t), dbH{t')\ , (53) 



The quality factor of the magnetization dynamics of 
most fcrromagnets is high (a ^ 0.01). Damping can 
therefore often be treated as a small perturbation. In 
the present Section we demonstrate that the damping ob- 
tained from linear response (Kubo) theory agrees^^ with 
that of the scattering theory of magnetization dissipation 
in this limit. At sufficiently low temperatures or strong 
elastic disorder scattering the coupling to phonons may 
be disregarded and is not discussed here. 

The energy dissipation can be written as 



E 




(48) 



where () denotes the expectation value for the non- 
equilibrium state. We are interested in the adiabatic 
response of the system to a time-dependent perturba- 
tion. In the adiabatic (slow) regime, we can at any time 
expand the Hamiltonian around a static configuration at 
the reference time t = 0, 




(49) 



i(r)-i>mBt(r) 



The static part, Hst, is the Hamiltonian for a magneti- 
zation for a fixed and arbitrary initial texture nist, as, 
without loss of generality, described by the collective 
coordinates S,a- Since we assume that the variation of 
the magnetization in time is small, a linear expansion in 
terms of the small deviations of the collective coordinate 
S£,i{t) is valid for sufficiently short time intervals. We can 
then employ the Kubo formalism and express the energy 
dissipation as 



where is the expectation value for the wave functions 
of the static configuration. Focussing on slow modula- 
tions we can further simplify the expression by expand- 



ing 



sut')~sut) + it' -t) 6Ut), 



(54) 



so that 



daH) = (daH) + dt'xab{t-t')S^b{t)+ 



dt'xab{t-t') {t' ~t) SUt)- 



(55) 



The first two terms in this expression, (daH) st + 
J_^dt'xab{t — t')S£_b{t), correspond to the energy vari- 
ation with respect to a change in the static magnetiza- 
tion. These terms do not contribute to the dissipation 
since the magnetic excitations are transverse, rh • m = 0. 
Only the last term in Eq. (|55p gives rise to dissipation. 
Hence, the energy loss reduces to^ 



E = i 



6iaS(b 



(56) 



uj=0 



where x^abi^) = JZ,dt [xab{t) + Xbait)] e^^* /2. The 
symmetrized susceptibility can be expanded as 



dH 



(50) 



i(r)->m3t(r) 



where the expectation value of the out-of-equilibrium 
conservative force 




daH 



(51) 



m(r)-!-mst(r) 



(/n - /m) {n\daH\m){m\dbH\n) + {a ^ h) 
2 hw + ii]- (e„ - Cm) 



(57) 

where |n) is an eigenstate of the Hamiltonian Hgt with 
eigenvalue e„, /„ = /(cn), /(e) is the Fermi-Dirac distri- 
bution function at energy e, and is a positive infinites- 
imal constant. Therefore, 
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(58) 



and the dissipation tensor 

de 



Tab = T^Y.[ -ir) {AdaH\m){m\dbHW)5{en ~ e™). (59) 



We now demonstrate that the dissipation tensor obtained 
from the Kubo Unear response formula, Eq. ([5^]) . is 
identical to the expression from scattering theory, Eq. 
PO)) . following the Fisher and Lee proof of the equiv- 
alence of linear response and scattering theory for the 
conductance<^ 

The static Hamiltonian HstH) = Hq + V{^) can be 
decomposed into a free-electron part Hq = —fi^W^ /2m 
and a scattering potential V{£,). The eigenstates of Hq 
are denoted |(/5s,q(e)) , with eigenenergies e, where s = ± 
denotes the longitudinal propagation direction along the 
system (say, to the left or to the right), and q a trans- 
verse quantum number determined by the lateral con- 
finement. The potential F(^) scatters the particles be- 



tween the propagating states forward or backward. The 
outgoing {+) and incoming (— ) scattering eigenstates 

of the static Hamiltonian Hgt are written as 'ipi^q\e)^, 

which form another complete basis with orthogonality re- 



lations ( ipi^"^ (e) 



d,^,-S,,,.S{e-e')M These 



wave functions can be expressed as 

gI^^V] \ips,q), where the retarded (-I-) and advanced (— ) 
Green's functions read Glf\e) = {e ± irj — H^t)^^- By 
expanding Tat in the basis of outgoing wave functions, 
l^tg) , the energy dissipation ([59| becomes 



Tab = TT ^ 

sq,s' q' 



de 



de 



daH 



(60) 



where wave functions should be evaluated at the energy e. 

Let us now compare this result, Eq. (j60p . to the direct scattering matrix expression for the energy dissipation, 
Eq. (|40p . The S'-matrix operator can be written in terms of the T-matrix as >S'(e;^) = 1 — 27rir(e;^), where the 
T- matrix is defined recursively by T = y[l + Gg^'^T]. We then find 



dT 



1 + VGit^ 



daH 



1 + Git^V 



The change in the scattering matrix appearing in Eq. (I4QD is then 



s' q' .sq 



-2ni {ifis^ql 



1 + VGit^ 



daH 



1 + Git^V 



daH 



(61) 



Since 



VI. CONCLUSIONS 



^sq.s'q' 

s'q' 



V'l^(£) 



(62) 



and SS'^ = 1, we can write the linear response result, 
Eq. ([60)) . as energy pumping (|40p . This completes our 
proof of the equivalence between adiabatic energy pump- 
ing in terms of the S'-matrix and the Kubo linear response 
theory. 



We have shown that most aspects of magnetization 
dynamics in ferromagnets can be understood in terms of 
the boundary conditions to normal metal contacts, i.e. 
a scattering matrix. By using the established numerical 
methods to compute electron transport based on scatter- 
ing theory, this opens the way to compute dissipation in 
ferromagnets from first-principles. In particular, our for- 
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malism should work well for systems with strong elastic 
scattering due to a high density of large impurity poten- 
tials or in disordered alloys, including Nii-ajFe^; (x ~ 0.2 
represents the technologically important "permalloy" ) . 

The dimensionlcss Gilbert damping tensors (j4ip for 
macrospin excitations, which can be measured directly 
in terms of the broadening of the ferromagnetic reso- 
nance, have been evaluated for Nii_a;Fea; alloys by ah ini- 
tio methods4^ Permalloy is substitutionally disordered 
and damping is dominated by the spin-orbit interaction 
in combination with disorder scattering. Without ad- 
justable parameters good agreement has been obtained 
with the available low temperature experimental data, 
which is a strong indication of the practical value of our 
approach. 

In clean samples and at high temperatures, the 
electron-phonon scattering importantly affects damping. 
Phonons arc not explicitly included here, but the scat- 
tering theory of Gilbert damping can still be used for 
a frozen configuration of thermally displaced atoms, ne- 
glecting the inelastic aspect of scatteringJ^ 

While the energy pumping by scattering theory has 
been applied to described magnetization dampingji^ it 
can be used to compute other dissipation phenomena. 
This has recently been demonstrated for the case of 
current-induced mechanical forces and damping;^ with 
a formalism analogous to that for current-induced mag- 
netization torques fi^iii^ 
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Appendix A: Adiabatic Pumping 

Adiabatic pumping is the current response to a time- 
dependent scattering potential to first order in the time- 
variation or "pumping" frequency when all reservoirs are 
at the same electro-chemical potentiali^ A compact for- 
mulation of the pumping charge current in terms of the 
instantaneous scattering matrix was derived in Ref. l39l . 
In the same spirit, the energy current pumped out of the 
scattering re gion has been formulated (at zero tempera- 
ture) in Ref. [23. Some time ago, we extended the charge 
pumping concept to include the spin degree of free- 
dom and ascertained its importance in magnetoelectronic 
circuitsji^ More recently, we demonstrated that the en- 
ergy emitted by a fcrromagnct with time-dependent mag- 
netizations into adjacent conductors is not only caused 



by interface spin pumping, but also reflects the energy 
loss by spin-flip processes inside the ferromagnet^ and 
therefore Gilbert damping. Here we derive the energy 
pumping expressions at finite temperatures, thereby gen- 
eralizing the zero temperature results derived in Ref. [l^ 
and used in Ref. [TBI . Our results differ from an earlier ex- 
tension to finite temperature derived in Ref. and we 
point out the origin of the discrepancies. The magneti- 
zation dynamics must satisfy the fluctuation-dissipation 
theorem, which is indeed the case in our formulation. 

We proceed by deriving the charge, spin, and energy 
currents in terms of the time dependence of the scattering 
matrix of a two-terminal device. The transport direction 
is X and the transverse coordinates are q = {y,z). An 
arbitrary single-particle Hamiltonian can be decomposed 
as 

Hir)^-^^+H^{x,e), (Al) 
where the transverse part is 

V{q) is an elastic scattering potential in 2 x 2 Pauli 
spin space that includes the lattice, impurity, and 
self-consistent exchange-correlation potentials, including 
spin-orbit interaction and magnetic disorder. The scat- 
tering region is attached to perfect non-magnetic electron 
wave guides (left a — L and right a — R) with constant 
potential and without spin-orbit interaction. In lead a, 
the transverse part of the 2x2 spinor wave function 
Lp"a\x,Q) and its corresponding transverse energy ei"' 
obey the Schrodinger equation 

ff±(e)^L")(e) = e(rVL"^(e), (A3) 

where n is the spin and orbit quantum number. These 
transverse wave guide modes form the basis for the ex- 
pansion of the time-dependent scattering states in lead 
a ~ L,R: 

Vv5L"^(e)e'^^'=^e-""""*/''ci"'='^\ (A4) 



where Cq annihilates an electron in mode n incident 
[a = 4-) or outgoing [a ~ ~) in lead a. The field opera- 
tors satisfy the anticommutation relation 

The total energy is ei"'^^ = h?k'^/2m + ei"-*. In the leads 
the particle, spins, and energy currents in the transport 
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direction are 

^ = ^ / ^^Tr, U-T- — ^ 



fie) 



h f ^ ^ / ft f 

— / dgTr, ^-tcr- —(T'f 

Arm / \ oa; aa; ; 



(A5a) 
(A5b) 



Ami 



9a; dx 



+H.C. 



(A5c) 



where we suppressed the time t and lead index a, cr = 
(cx, Cy, cTz) is a vector of Pauh matrices, and Tr^ denotes 
the trace in spin space. Note that the spin current Is 
flows in the x-direction with polarization vector Is/Is- 
To avoid dependence on an arbitrary global potential 
shift, it is convenient to work with heat /^'^^ rather than 
energy currents J^'^-' : 

i(i^t) ^ i'-'^t) - fii'^p\t) , (A6) 

where fi is the chemical potential. Inserting the waveg- 



uide representation (|A4p into (jASp . the particle current 
reads^ 



Attiti 



POO 

/ dkdk' V {ak + a'k') 
Jo I 



We are interested in the low- frequency limit of the Fourier 
transforms /^''^(w) = j^^dte"^Hl^\t). Following Rcf. 
I41I we assume long wavelengths such that only the inter- 
vals with k k' and a = a' contribute. In the adiabatic 
limit w ^ this approach is correct to leading order in 
hjj/ep, where is the Fermi energy. By introducing the 
(current-normalized) operator 



^{n<7)i {nk)\ _ 



1 



~(nka) 
- ^1-^ 



(A8) 



dk 



which obey the anticommutation relations 



The charge current can be written as 

1 1"°° 



(A9) 



(AlO) 



We operate in the linear response regime in which applied 
voltages and temperature differences as well as the exter- 
nally induced dynamics disturb the system only weakly. 
Transport is then governed by states close to the Fermi 
energy. We may therefore extend the limits of the en- 
ergy integration in Eq. (jAlOp from (ei"'',oo) to (—00 
to 00). We relabel the annihilation operators so that 



.(nfc) 



^(nk) 



denotes particles incident on the scattering 

region from lead a and c^"^' denotes particles 

leaving the scattering region by lead a. Using the Fourier 
transforms 



— 00 

1 

2t^ 



decS^"'\e)e-"'*/^, 



we obtain in the low-frequency limitsi 



(All) 
(A12) 

(A13) 



r 



all wave-guide modes {S^"''}. Analogous calculations lead 
to the spin current 



and the energy current 



,i dda dba 



H.C.. 



(A14) 



(A15) 



Next, we express the outgoing operators b{t) in terms 
of the incoming operators a{t) via the time-dependent 
scattering matrix (in the space spanned by all waveguide 
modes, including spin and orbit quantum number): 

= E / dt'S^p{t,t')dp{t'). (A16) 



where ba is a column vector of the creation operators for 



When the scattering region is stationary, Sai3{t,t') only 
depends on the relative time difference t — t', and its 
Fourier transform with respect to the relative time is 
energy independent, i.e. transport is elastic and can 
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be computed for each energy separately. For time- 
dependent problems, Sapit^t') also depends on the total 
time t + t' and there is an inelastic contribution to trans- 
port as well. An electron can originate from a lead with 
energy e, pick up energy in the scattering region and end 
up in the same or the other lead with different energy e'. 

The reservoirs arc in equilibrium with controlled lo- 
cal chemical potentials and temperatures. Wc insert the 
S'-matrix (|A16|) into the expressions for the currents, 



Eqs. (|A13|) . (|A14|) . (|A15|) . and use the expectation value 
at thermal equilibrium 

(4(")(i2)al;"^(tl)) ^Sn,nS^pUh^t2)/27Th, (Al7) 
\ / cq 

where ffs{ti-t2) = i2Trh)-^ J de-^<*^-*^'>/^ fa{e) and 
/c((e) is the Fermi-Dirac distribution of electrons with 
energy e in the a-th reservoir. Wc then find 



2TTh(bi{t)(Tb^{t)'j^^^J2 J dtldt2S*^p{t,t2)eTSap{t,h)fp{ti-t2), 



(A18) 
(A19) 
(A20) 



Next, we use the Wigner representation (jBl 



S{t,t') 



2Trh 



deS 



t + t' 



(A21) 



and by Taylor expanding the Wigner represented S-matrix S{{t + t')/2,e) around S{t,e), S{{t + t')/2,e) 
Er=o d^S{t, e)it' ~ i)'V(2""!), wc find 



ZTTAl 



(A22) 



and 



hdtS{t, t') = r dee-'^(*-*')/^e*'*^=^'/2 (\^q^ _ ^\ g^^^ ^^23) 
2nhJ_^ \2 J 

The factor 1/2 scaling the term hdtS{t, e) arises from commuting e with e''''^^^*/^. The currents can now be evaluated 



(A24a) 
(A24b) 



1 

1 



de 



(A24c) 



where the adjoint of the S-matrix has elements S^^^ = 5*^"'" ■* . 

We are interested in the average (DC) currents, where simplified expressions can be found by partial integration 
over energy and time intervals. We will consider the total DC currents out of the scattering region, /(°"') = — /„, 
when the electrochemical potentials in the reservoirs are equal, /^(e) = /(e) for all a. The averaged pumped spin and 
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energy currents out of the system in a time interval r can be written compactly as 



r(c) _ 







dt / deTr 



1 



t(s) 

™' 27rfir Jo 



1 



27rftr 



1 



2-KhT 



dt J deTi <, <j 
dt I deTi- 
dt / deTr 



ih d\ 
ih d 



- ef{e) 



2 dt 



(A25a) 
(A25b) 

(A25c) 



where Tr is the trace over all waveguide modes (spin 
and orbital quantum numbers). As shown in Ap- 
pendix [C] the charge pumped into the reservoirs vanishes 
for a scattering matrix with a periodic time dependence 
when, integrated over one cycle: 

4Si=0. (A26) 

This reflects particle conservation; the number of elec- 
trons cannot build up in the scattering region for peri- 
odic variations of the system. We can show that a similar 
contribution to the energy current, i.e. the first line in 
Eq. (|A25cp . vanishes, leading to to the simple expression 

(A27) 

Expanded to lowest order in the pumping frequency the 
pumped spin current (jA25bp becomes 

(A28) 

This formula is not the most convenient form to com- 
pute the current to specified order. SS"^ also contains 
contributions that are linear and quadratic in the pre- 
cession frequency since S{t, e) is the 5'-matrix for a time- 
dependent problem. Instead, we would like to express the 
current in terms of the frozen scattering matrix Sh{t, e). 
The latter is computed for an instantaneous, static elec- 
tronic potential. In our case this is determined by a mag- 
netization configuration that depends parametrically on 
time: S'fi.(i, e) = 5[m(t),e]. Using unitarity of the time- 



dependent ^-matrix (as elaborated in Appendix ICj) . ex- 
pand it to lowest order in the pumping frequency, and 
insert it into (jA28|) leads to^S 

(A29) 

We evaluate the energy pumping by expanding (|A27[) 
to second order in the pumping frequency: 

Tie) h r dt [ ( dS^ IdSdS^X 

(A30) 

As a consequence of unitarity of the S'-matrix (see Ap- 
pendix [C]), the first term vanishes to second order in the 
precession frequency: 

where, at this point, we may insert the frozen scattering 
matrix since the current expression is already propor- 
tional to the square of the pumping frequency. Further- 
more, since there is no net pumped charge current in 
one cycle (and we are assuming reservoirs in a common 
equilibrium) , the pumped heat current is identical to the 

pumped energy current, = 1^"^^. 

Our expression for the pumped energy current (|A31[) 
agrees with that derived in Ref. 1271 at zero temperature. 
Our result (|A3ip differs from Ref. |40 at finite tempera- 
tures. The discrepancy can be explained as follows. In- 
tegration by parts over time t in Eq. (|A27p , using 



ih d 
'2dt 



ih 



'~dt 



5t 



ef 



and the unitarity condition from Appendix [Cl 



ih d 
'2dt 



S'' -2 



'lide 



ih d\ f ih d 

^^Ydir' Tat 



s 



5t = 



ih d 
Ydt 



ih d 

'2 at 



s 



si 



(A32) 



(A33) 
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the DC pumped energy current can be rewritten as 



1 



dt 



&Tr 



efie) 



Next, we expand this to the second order in the pumping frequency and find 

4^1 = ^ r? /*Tr(./(e)(55t-l)-<5J)!^^5t_,(a^./)!^g^t 



2 dt 



8 dt^ 



(A34) 



(A35) 



r 



This form of the pumped energy current, Eq. (jA35p . 
agrees with Eq. (10) in Ref. HO if one {incorrectly) as- 
sumes SS^ = 1. Although for the frozen scattering ma- 
trix, ShSf^ = 1, unitarity does not hold for the Wigncr 
representation of the scattering matrix to the second or- 
der in the pumping frequency. (5" 5*^ — 1) therefore does 
not vanish but contributes to leading order in the fre- 
quency to the pumped current, which may not be ne- 
glected at finite temperatures. Only when this term is 
included our new result Eq. (|A31|) is recovered. 



Appendix B: Fourier transform and Wigner 
representation 

There is a long tradition in quantum theory to trans- 
form the two-time dependence of two-operator correla- 
tion functions such as scattering matrices by a mixed 
(Wigner) representation consisting of a Fourier transform 
over the time difference and an average time, which has 
distinct advantages when the scattering potential varies 
slowly in time4i In order to establish conventions and no- 
tations, we present here a short exposure how this works 
in our case. 

The Fourier transform of the time dependent annihi- 
lation operators are defined in Eqs. (jAll[) and (jA12l) . 



Consider a function A that depends on two times ti 
and t2, A = A{ti,t2)- The Wigner representation with 
t ~ [ti + t2)/2 and t' = ti — t2 is defined as: 



AitiM) = A / deA{t,e) 
A(t,e) = 



-ie{ti—t2)/h 



(Bl) 



t' t' 



dt'A[t+^-,t-^-]e'''''^. (B2) 
We also need the Wigner representation of convolutions, 

/oo 
dt'A{ti,t')B{t',t2). (B3) 
-oo 

By a series expansion, this can be expressed as^ 

{A (E) B){t, e) = e-<^^^^-^^''^^)^^A{t, e)B{t, e) (B4) 

which we use in the following section. 

Appendix C: Properties of S-matrix 

Here we discuss some general properties of the two- 
point time-dependent scattering matrix. Current conser- 
vation is reflected by the unitarity of the S'-matrix which 
can be expressed as 



(3n' s' 



(Cl) 



Physically, this means that a particle entering the scattering region from a lead a at some time t is bound to exit the 
scattering region in some lead /3 at another (later) time t' . Using Wigner representation (jBl|) and integrating over 
the local time variable, this implies (using Eq. (jB4[) ') 



1 



{S®S^) (t,e) = e-<^'^''-^''^')/'5(t,e)5t(i^g)^ (C2) 



where 1 is a unit matrix in the space spanned by the wave guide modes (labelled by spin s and orbital quantum 
number n). Similary, we find 

1={S^® S) {t,e) = e+*(^"^"'-^"'^")/'5t(t,e)5(t,e). (C3) 
To second order in the precession frequency, by respectively subtracting and summing Eqs. (|C2[) and (|C3p give 
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and 



Tr{S'5^-l} =Tr 



dt^ 9e2 ~ dtde dtde ^ 9e2 df^ 



(C5) 



Furthermore, for any energy dependent function Z[e) and arbitrary matrix in the space spanned by spin and transverse 
waveguide modes Y , Eq. (|C2I) implies 



dt / deZ{e)Tr 



(C6) 



Integration by parts with respect to t and e gives 



dt / deTr 



(C7) 



which can be simphfied to 



dt / deTr 



Similarly from (jC3p . we find 



(C8) 



IJ^dtJ deTrl^(^SHt,e 



-1]Y'>=0 



(C9) 



r 



Using this result for y = 1 and Z{e) = /(e) in the 
expression for the DC particle current (|A25al) . we see 
that unitarity indeed implies particle current conserva- 
tion, la"^ = for a time-periodic potential. However, 
such a relation does not hold for spins. Choosing Y = cr, 
we cannot rewrite Eq. (|C9p in the form (|A25bp . unless 
the iS-matrix and the Pauli matrices commute. Unless 
the 5- matrix is time or spin independent, a net spin cur- 
rent can be pumped out of the system, simultaneously 
exerting a torque on the scattering region. 



Furthermore, choosing Z{() = de ' /{e'), Y ~ 1 and 
expanding the difference between (|C9I) and (jCSp to sec- 
ond order in frequency gives 

^J\tJdeTr!^fie)^^S\t,e)^^0, 

which we use to eliminate the first term in the expression 
for the energy pumping, Eq. (|A30p . 
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